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摘 X8: 蒸 散发 (ET) 是 区 域 能 量 平衡 以 及 水 量 平衡 的 关键 环节 ， 精 确 估 算 蒸 散发 ， 对 于 提高 水 分 利用 效率 以 
及 优化 区 域 用 水 结构 具有 重要 意义 ， 而 冠 层 阻抗 则 是 准确 估算 蒸 散发 的 一 个 重要 变量 。 为 了 确定 冠 层 阻 抗 模 
型 区 域 适 用 性 、 解 决 其 参数 化 问题 ， 本 研究 基于 黑河 重大 研究 计划 已 有 的 通 量 观测 数据 ， 以 Irmak 模型 为 基础 ， 
考虑 微 气象 因子 与 冠 层 阻抗 之 间 的 关系 ,增加 了 大 气 CO, 浓度 对 冠 层 阻抗 的 影响 ,构建 了 未 考虑 COS 和 考虑 
CO; 影响 的 两 种 Irmak 模型 ， 并 将 其 与 Penman-Monteith(P-MD) 模 型 耦合 ， 利 用 已 有 涡 度 相关 数据 ， 分 析 和 检验 
了 两 种 冠 层 阻抗 模型 对 环境 变量 和 大 气 CO; 浓度 响应 的 模拟 结果 ， 并 对 模型 参数 进行 敏感 性 分 析 。 结果 表明 : 
将 考虑 大 气 CO, 浓度 影响 的 Irmak 模型 与 Penman-Monteith 模型 看 合 ， 能 够 更 好 地 模拟 玉米 冠 层 阻抗 和 蒸 散 量 
对 外 部 环境 变量 的 响应 过 程 。 在 参数 率 定 期 该 模型 所 模拟 的 冠 层 阻抗 和 蒸 散 量 与 实测 值 之 间 的 RO 分别 达 0.76 
fu 0.95, RMSE 分 别 达 33.1 sm 和 34.5 W-m ^; 模型 验证 期 冠 层 阻 抗 和 蒸 散 量 模拟 值 与 实测 值 之 间 的 R? 分 别 
3k 0.68 和 0.90, RMSE 分 别 达 63.2 sm ! fu 49.0 W:m?。 两 个 独立 验证 点 结果 表明 考虑 了 大 气 CO, 浓度 影响 的 
Irmak 模型 具有 较 好 的 空间 可 移植 性 和 适应 性 ， 模 型 能 够 较为 准确 地 模拟 玉米 在 整个 生长 季 半 小 时 时 间 尺 度 
上 的 农田 耗 水 过 程 。 敏 感性 分 析 表 明 玉 米 冠 层 阻抗 及 其 蒙 散 量 对 净 辐 射 和 相对 湿度 变化 最 为 敏感 ， 其 次 是 气 
温 、 叶 面积 指数 和 大 气 CO; 浓度 。 本 文 所 构建 的 考虑 大 气 CO 浓度 对 于 玉米 冠 层 阻抗 影响 的 Irmak 模型 能 够 
较为 准确 地 估算 作物 蒸 散 量 ， 并 可 为 种 植 结 构 调 整 、 土 地 利用 方式 改变 以 及 大 气 CO 浓度 变化 环境 下 的 农田 
耗 水 研究 提供 一 定 的 研究 依据 。 
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Abstract: Evapotranspiration (ET) is critical for energy and water balance in agricultural systems. Accurate estimation or 
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measurement of ET is therefore important in improving water use efficiency and optimizing the structure of regional water use. 
Canopy resistance is one of the most important variables in the estimation of ET. The accuracy of simulation of the response process 
of canopy resistance to environmental variables is critically important for crop ET research. A convenient approach to simulate the 
response process of canopy resistance to multiple factors is based on the relationship between measured latent heat, climatic variables 
and by using the modified Penman-Monteith (P-M) equation. However, this method has certain limitations in some practical 
applications due to the lack of a few effective parameters. Another approach is to construct empirical and semi-empirical models 
using multiple factors (such as the Irmak model) based on measured data combined with the rearranged P-M equation. Based on 
canopy resistance values calculated by the rearranged P-M equation and on maize data (for the period May to September 2012) 
collected from the three eddy covariance observation stations in Heihe River Basin, this study constructed Irmak model taking into 
account the effect of atmospheric CO; concentration of half-hourly and daily time-steps to simulate the response processes of maize 
to environmental variables such as net radiation (R,), air temperature (T,), leaf area index (LAT), relative humidity (RH), wind speed 
(U3), aerodynamic resistance (ra), effective soil water content (0) and atmospheric CO» concentration. In the study, the performance of 
the two Irmak models were tested with measured values of latent heat from the eddy covariance systems of the other two verification 
points. Besides, the sensitivity of environmental variables was analyzed. The results indicated that the improved Irmak model which 
took into account the effect of atmospheric CO; concentration well estimated canopy resistance and ET. The coefficients of 
determination (R?) for canopy resistance and ET were respectively 0.76 and 0.95 for the calibration phase, with root mean square 
errors (RMSE) of 33.1 s:m™' and 34.5 W:m ?. Meanwhile, R? for canopy resistance and ET were respectively 0.68 and 0.90 for the 
validation phase, with RMSE of 63.2 sm ! and 49.0 W-m ?. The two verification points showed that the improved Irmak model had a 
good performance and strong regional applicability and spatial portability. The model also simulated the response processes of canopy 
resistance to environmental variables and reflected the effect of the variations in atmospheric CO» concentration on ET. Sensitivity 
analysis of the improved Irmak model showed that canopy and ET were the most sensitive to net radiation and relative humidity, 
followed by air temperature, leaf area index and atmospheric CO; concentration. The improved Irmak model used in this study was 
applicable in estimating crop water consumption and the accuracy of ET of maize, in providing scientific basis for improvements in 
water use efficiency and in optimizing the structure of regional water use under increased future atmospheric CO» concentration. 

Keywords: Canopy resistance; Evapotranspiration (ET); Penman-Monteith (P-M) model; Irmak model; Atmospheric CO; 


concentration 


植物 叶片 上 的 气孔 是 植物 体 进行 体内 外 气体 交 
换 的 通道 , 气孔 的 开 闭 行 为 会 影响 植物 体 的 蒸腾 以 
及 光合 作用 。 植 物 组 织 中 的 水 分 通过 蒸腾 作用 进入 
大 气 ， 从 而 实现 植被 与 外 界 的 能 量 交换 门 ,该 过 程 
受到 外 部 环境 因子 和 植物 生理 等 因素 的 共同 控制 和 
调节 个。 冠 层 阻抗 是 水 分 在 作物 冠 层 与 大 气 进行 循 
环 所 需要 克服 的 阻力 ,如 何 更 为 准确 地 模拟 冠 层 阻 
抗 ， 对 于 精准 估算 蒸 散发 ECT)， 进 而 对 于 提高 水 分 
利用 效率 和 优化 区 域 用 水 结构 县 有 重要 意义 1。 

许多 经 验 和 半 经 验 的 模型 被 用 来 估算 冠 层 阻抗 
对 环境 变量 的 响应 过 程 , 例如 Jarvis 模型 、 双 源 耦 合 
模型 、Irmak 模型 、S-W 模型 、K-P RESTA, 
其 中 最 为 典型 的 Jarvis 模型 将 冠 层 阻抗 看 成 是 外 部 
环境 变量 胁迫 的 结果 , 建立 了 形式 简单 灵活 的 冠 层 
阻抗 模型 "但 生理 意义 不 明确 ,未 考虑 各 因子 之 
间 的 相互 作用 , 模型 的 复杂 性 随 着 外 部 环境 变量 的 
增多 而 增加 ,模型 精度 随 着 模拟 时 间 的 延长 而 降低 ， 
尤其 是 最 小 冠 层 阻力 的 确定 较为 困难 "Li 等 
基于 石 羊 河流 域 大 量 试验 所 建立 的 耦合 了 土壤 蒸发 
和 植被 蒸腾 的 耦合 模型 ， 对 玉米 (Zea mays) 和 酿酒 
葡萄 (Fiis vinifera) 均 具有 较 好 的 模拟 效果 ,但 模型 


参数 较 多 , 个 别 参 数 不 易 获取 中 ， 部 分 参数 需要 试 
验 数据 进一步 修正 ， 模 型 精度 在 作物 叶 面积 指数 
(LAD 较 低 情况 下 仍 受 到 限制 。Irmak 等 所 利用 多 年 
观测 资料 , 考虑 微 气象 因子 与 冠 层 阻抗 之 间 的 关系 ， 
采用 广义 非 线 性 回归 方法 设计 了 7 个 模型 ,模型 的 
复杂 性 随 着 微 气象 因子 的 增多 而 增加 ， 该 模型 虽 为 
生理 意义 不 明确 的 经 验 模型 ， 但 形式 简洁 ， 计 算 简 
单 ， 使 用 方便 快捷 。Irmak 等 与 利用 其 构建 的 阻抗 模 
型 对 玉米 的 冠 层 阻抗 及 蒸 散 量 进行 模拟 ， 在 后 来 的 
研究 中 又 利用 该 模型 对 大 豆 (Glycine max) 的 冠 层 阻 
HARRE HITR, 模拟 精度 很 高 ， 然 而 模型 
的 适应 性 和 可 移植 性 有 待 于 进一步 研究 。 

众多 研究 成 果 呈 "7 表明 以 上 模型 均 具 有 较 好 
的 模拟 效果 以 及 较 高 的 适用 性 , 许多 新 模型 在 此 基 
础 上 进行 了 修正 和 改进 ,考虑 的 因子 也 更 加 全 面 , 但 
这 些 模型 忽略 了 大 气 CO, 浓 度 对 冠 层 阻抗 的 影响 1。 
现 有 气候 变化 条 件 下 ， 由 于 大 气 CO 浓度 日 变化 幅 
度 较 小 ， 对 冠 层 阻 抗 的 影响 轻微 , 但 从 长 期 模型 预 
测 的 角度 考虑 ， 尤 其 是 在 未 来 气候 和 生态 系统 水 磋 
平衡 变化 的 模拟 预测 中 , 考虑 大 气 CO, 浓度 变化 对 
蒸 散 量 的 影响 更 加 重要 。 现 有 试验 表明 大 气 CO 
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度 的 改变 会 影响 植物 体 与 外 界 水 气 通 量 交 换 ， 大 气 
CO; 浓度 增加 , 会 导致 许多 植物 的 气孔 导 度 下 降 王 ”1 
冠 层 阻抗 增加 ， 蒸 腾 作 用 减弱 ， 尤 其 是 高 浓度 CO; 
能 诱导 和 气孔 关闭 己 20。Wand 等 上 研究 发 现 CO; 浓 
EDERT, Cs 和 Cs 作物 的 气孔 导 度 分 别 降 低 29% 和 
24%。 在 全 球 变化 背景 下 ， 大 气 CO, 浓度 上 升 和 水 资 
源 短缺 已 经 对 区 域 农 业 的 可 持续 发 展 产生 了 重要 影 
ng pip IE CO; 浓度 对 冠 层 阻抗 的 影响 ， 对 于 
揭示 植物 光合 作用 和 耗 水 规律 ,提高 水 分 利用 效率 
BAEAN, 

本 研究 以 Irmak 模 型 为 基础 ， 基 于 Morison $0 
和 Easterling EPX F CO, 浓度 升 高 情景 下 作物 的 
响应 规律 研究 ， 假 设 大 气 CO 浓度 对 玉米 叶片 气孔 
导 度 的 影响 呈 线 性 趋势 ， 考虑 叶 面积 指数 (LAT)、 净 
辐射 (R,)、 相 对 湿度 (RH)、 气 温 (7,)、 土 壤 有 效仿 水 
量 (0)、 空 气动 力学 阻抗 (,)、3 m 高 度 处 风速 (U3) 和 
XA CO, 浓度 等 与 冠 层 阻抗 之 间 的 关系 ,构建 了 考 
虑 CO; 浓度 和 未 考虑 CO, 浓 度 影响 的 冠 层 阻抗 模型 ， 
并 将 其 与 Penman-Monteith 模型 耦合 ， 确 定 模型 在 
区 域 上 的 适用 性 及 其 参数 化 问题 。 利 用 “黑河 流域 生 
态 -水 文 过 程 集成 研究 ”重大 研究 计划 已 有 的 通 量 观 
测 数 据 ， 建 立 冠 层 阻 抗 与 环境 因子 之 间 的 关系 ， 重 
点 考虑 了 大 气 CO 浓度 对 冠 层 阻 抗 的 影响 ， 筛 选 出 


能 够 反映 CO; 浓度 变化 环境 下 的 冠 层 阻 抗 模型 ， 以 
期 为 揭示 作物 耗 水 规律 ,精准 估算 作物 蒸 散 量 ， 提 
高 水 分 利用 效率 以 及 作物 蒸 散 量 在 未 来 大 气 CO; 浓 
度 上 升 条 件 下 的 响应 机 制 提供 一 定 的 科学 依据 。 


1 研究 区 概况 、 数 据 处 理 与 模型 简介 
11 ”研究 区 概况 

黑河 流域 (图 1) 位 于 河西 走廊 中 部 ， 地 处 青藏 高 
原 与 内 蒙古 高 原 的 过 渡 地 带 ， 是 我 国 西北 干旱 半 干 
旱 区 第 二 大 内 陆 河流 域 ， 流 域 面积 约 04.3x10* km 
该 流域 气候 和 干燥， 降水 稀少 , 水 资源 分 配 不 均 。 流 
域 主 要 分 为 上 游 和 姥 连 山地 ， 中 游 走 廊 绿 洲 平原 以 及 
下 游 阿 拉 善 高 原 3 个 地 貌 类 型 区 。 黑 河中 游 是 指 黑 
河 干流 出 山口 葛 落 峡 以 下 至 正义 峡 之 间 的 地 势 平 
坦 区 域 , 为 典型 的 温带 大 陆 性 气候 , 年 降水 量 为 
116.8 mm, FREE 2365.6 mm, 年 平均 气温 
7.6 Co 7 月 到 9 月 的 降水 量 占 全 年 降水 量 的 60% 
左右 ,冬季 降水 量 仅 占 3%, 年 降水 量 随 着 海拔 高 
度 增加 和 距 河 流 距离 增加 而 逐渐 减少 。 年 日 照 时 数 


3085 h, 无 霜 期 165 d。 中 游 绿 洲 区 耕地 面积 占 全 
流域 耕地 总 面积 的 95%， 水 资源 消耗 占 全 流域 的 
68%， 农 田 蒸 散 为 水 资源 消耗 的 主要 项 ， 玉 米 为 当 
地 主要 农作物 之 一 站。 


| EI 玉米 Maize 1 


站 点 EC station 


积 雪 Snow cover — 黑河 Heihe River 
= 首 蒂 Alfalfa [73 建设 用 地 Construction land 
国 草地 Grassland B thi Other farm land 


B 阔 叶 林 Broadleaved forest | KÆ Barley 

= | HEB Bare land Bl 水 体 Water body 
[93 针 叶 林 Coniferous forest B 混 地 wet land 
棉花 Cotton [33 灌 从 地 Bush fallow 


流域 试验 站 点 及 其 土地 利用 


1 GER 
Fig. 1 Experimental sites and its’ land use types in Heihe River Basin 


http://www.ecoagri.ac.cn 


250 中 国生 态 农业 学 报 2017 58-2578 
1.2 ”数据 来 源 与 处 理 CO 

=1+|1- 一 2 |xA 4 

本 研究 所 采用 的 自动 气象 站 、 通 量 、 作 物 株 高 f(cox) 4 3 gLC0, e 


等 数据 均 来 源 于 黑河 计划 数据 管理 中 心 (http:// 
heihedata.org/)， 通 量 观测 和 矩阵 中 的 Site 8(100922'35"E, 
38°52'21"N; 1550.06 m) 用 于 模型 参数 的 率 定 ，Site 
11(100920'31"E, 38?52'12"N; 1575.65 m)、 大 满 站 
(100922'20"E, 38?51'20"N; 1556.06 m) 的 数据 用 于 
模型 的 验证 (图 1), 每 个 通 量 观测 系统 主要 由 三 维 超 
声 风 温 仪 (CSAT3，Campbell Scientific，USA) 和 开路 
COXH2O 红外 气体 分 析 仪 (Li-7500A, Li-Cor Inc., USA) 
组 成 , 原始 数据 采样 频率 为 10 Hz, 8 30 min 输出 1 
组 平均 通 量 值 以 及 10 Hz 原始 数据 "1。3 个 通 量 观测 
站 点 下 垫 面 作物 均 为 玉米 ， 玉 米 种 植 垄 距 为 50.8 cm, 
行距 为 43.3 cm, 株距 为 22 cmo. 玉米 株 高 由 实测 数 
据 进 行 3 次 样 条 插值 内 插 到 每 半 小 时 ; 叶 面 积 指数 
(LAD) 来 源 于 美国 国家 航天 局 (NASA) 发 布 的 MCD15A3H 
叶 面 积 指数 产品 (https://search.， earthdata.nasa.gov/), 
该 数据 为 4 d 合成 的 500 m DIRAY L4 级 叶 面积 指 
数 数 据 ,， 采 用 3 次 样 条 插值 法 均匀 内 插 到 每 天 ， 并 
假设 同一 天 不 同时 刻 叶 面 积 指数 不 发 生变 化 。 研 究 
时 间 序 列 为 整个 作物 生长 季 数 据 (2012 年 5 月 初 一 
2012 年 9 月 底 )， 时 间 段 为 9:00— 18:30。 原始 涡 度 相 
关 仪器 信息 及 数据 处 理 见 文献 [31]， 原 始 涡 度 相关 数 
据 已 经 过 野 点 剔除 、 延 迟 时 间 校 正 、 坐 标 旋 转 、 响 
应 频率 修正 和 严格 的 质量 控制 等 步骤 ,但 还 需要 对 
部 分 原始 数据 进行 再 处 理 55， 处 理 规则 如 下 : 1) 对 
原始 涡 度 相关 数据 进行 能 量 闭合 度 检查 ,将 能 量 闭 
合 度 超出 0.5~1.5 之 间 的 数据 剔除 ， 其 余 能 量 未 闭合 
的 数据 进行 强制 闭合 ,闭合 方法 见 [33]; 2) 去 除 有 空 
值 的 数据 ; 3) 剔 除 异 常数 据 ， 如 实测 的 潜 热 通 量 
(LE), re oo>2 000 或 mwr<00。 ow 为 利用 Penman-Monteith 
公式 反 推 的 冠 层 阻抗 值 , 此 处 定义 为 实测 值 )。 
13 ”模型 简介 
1.3.1 Irmak 模型 

基于 Irmak MHEARA eH 
(LAT)、 净 辐射 (R,)、 相 对 湿度 (RH)、 气 温 (7,)、 土 壤 
有 效 合 水 量 (0)、 空 气动 力学 阻抗 (x,)、3 m 高 度 处 风 
速 (3) 和 大 气 CO; 浓度 (CO,) 等 外 部 环境 变量 对 冠 层 
阻抗 的 影响 ,构建 了 有 Co 浓度 影响 和 无 CO, 浓度 
影响 的 冠 层 阻抗 模型 ， 形 式 如 下 : 
ri 7 exp[a * bR, *cT, - dRH * eU, + gr, * hLAI +if (0)] 


(1) 
a *bR, c cT, +dRH+ 
“2 = apl so, + gr, en (2) 


f(8) -(8-8,) 8; —8,) 3) 


AA: ron 和 rn 分别 为 在 Irmak 模型 的 基础 上 未 考 
EAS CO, 浓度 影响 和 考虑 了 大 气 CO 浓度 影响 的 
冠 层 阻抗 模型 ; /(9) 为 土壤 合 水 量 ; f(CO;)79 CO; 
浓度 胁迫 画 数 ， 式 (4) 中 Ag; CO; 为 CO, 浓度 增加 1 
fij, 叶片 气孔 导 度 减 小 的 倍数 路 ， 本 研究 取 0.3; 
0 为 田间 持 水 量 ， 取 0.34 cmm”; 0, HERRA, 
取 0.1 cm3.cm-3; a. b. cx d. ex g h i23979 8E 
定 的 经 验 系 数 。 
1.3.3 Penman-Monteith 模型 
Penman-Monteith 模型 是 Monteith 在 Penman 模 

型 的 基础 上 发 展 而 来 , 模型 如 下 : 

A(R, -G)* C, p, VPD/r, 


AET (5) 

Atytyxr.lr, 
由 公式 (5) 可 以 得 出 : 
C. xp, x VPD 
n| A(R, -G)«— A AET x (A +7) 
r, 
ji = a 6 
"c ob yxAET (9 


式 中 : x 为 利用 Penman-Monteith 公式 反 推 得 到 的 
冠 层 阻抗 ; 4 为 汽化 潜 热 , [kg …; ET 为 实际 蒸 散 量 ， 
Wm’; A 为 饱和 水 汽 压 与 温度 关系 曲线 斜率 , kPa C; 
RR, 为 净 辐 射 , Wm; GARAR, Wm’; CGA 
空气 定 压 比 热 , Jk C; p, 为 空气 密度 ，kgm 7; 
VPD 为 饱和 水 汽 压 差 , kPa; 7 为 冠 层 阻抗 , sm n 
为 空气 动力 学 阻抗 , sm '; y IPERE, kPa C, Æ 
气动 力学 阻抗 计算 如 下 m1: 

enm 
_ Zom Zoh (7) 

K^U, 
式 中 : x; 为 空气 动力 学 阻抗 , sm; Za ZA GS S, 
m; Z, 为 测 湿度 高 度 , m; d 为 震 平 面 位移 , m, d—2/3 h; 
有 为 作物 高 度 , m; Zu 为 控制 动量 传递 的 粗糙 长 度 , m, 
Zom=0.123 h; Zu 为 控制 热 通 量 和 水 汽 传输 的 粗糙 长 
E, m, Zon =0.1 Zom; K 为 卡 曼 常数 , 0.41; UV, H Z AE 
处 的 风速 , ms; 本 文中 Z,-3 m, Z;-5 mo 
14 ”模型 评价 
本 文采 用 决定 系数 (R*”), 均值 偏 移 误 差 (mean bias 
error,，MBE) 以 及 均 方 根 误差 (root mean square error, 
RMSE) 来 评价 模型 模拟 值 (6%) 与 实测 值 (O06。) 之 间 的 
差异 , 检验 模型 模拟 精度 。 
2 ducha 

aO 7 Oa * Gs O3] 


a 


R? =1 


(8) 
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] 4 
MBE = - 2 En ES Os) (9) 


i-l 


RMSE = Se -0x } (10) 
ni 


式 中 : Ea DRA, Os 73 323, Op 为 实测 值 的 
平均 值 , n 为 样本 容量 。 本 文中 约定 利用 涡 度 相关 实 
测 的 潜 热 通 量 并 结合 Penman-Monteith 公式 反 推 得 
到 的 冠 层 阻 抗 为 玉米 冠 层 阻抗 的 实测 值 (x, o), AA 
模型 拟 合 的 冠 层 阻 抗 为 模拟 值 (x. m 利用 涡 度 相关 
实测 的 潜 热 通 量 为 ET 的 实测 值 (ET o) 利用 模型 模 
拟 的 ET 为 模拟 值 (ET m)o 


2 ”模型 模拟 效果 与 验证 
21 ”模型 参数 率 定 及 验证 
本 研究 将 Site 8 点 的 数据 (数据 时 段 为 2012 年 6 


月 6 日 一 9 月 20 日 ) 分 成 前 后 两 部 分 ,前 一 部 分 数据 的 
观测 时 段 为 6 月 6 日 一 8 月 5 A, 用 于 模型 参数 率 定 ; 
后 一 部 分 数据 的 观测 时 段 为 8 月 6 日 一 9 H 20 A, 用 
于 模型 参数 验证 ,将 参数 率 定期 的 数据 分 别 代 入 所 构 
建 的 两 个 Irmak 模型 中 , 采用 最 小 二 乘法 进行 非 线 性 
回归 ， 率 定 模 型 参数 ， 建 立 冠 层 阻抗 与 外 部 环境 因 
子 之 间 的 经 验 公式 ; 并 采用 验证 期 的 数据 对 所 建立 的 
模型 进行 验证 , 最终 得 到 相关 参数 的 优选 值 见 表 1。 为 
了 进一步 检验 模型 在 整个 作物 生育 期 内 的 模拟 精度 ， 
采用 Site 11 和 大 满 两 个 站 点 作为 进一步 的 验证 点 。 
2.1.1 Irmak 模型 及 其 改进 形式 的 模拟 效果 

利用 所 构建 的 未 考虑 CO, 浓度 (x。 1) 和 考虑 
CO;(r, D) 的 两 个 模型 ， 分 别 对 Site 8 的 玉米 冠 层 阻 
抗 及 切 散 量 进行 模拟 , 模型 相关 参数 见 表 1， 模 型 对 
比 结果 见 图 2 和 图 3。 


表 1 RE CO, REO. rm 和 考虑 CO 浓度 (rs pz) 的 Irmak 模型 中 相关 参数 的 最 优 值 


Table 1 Optimum values of parameters of the two Irmak models non-considering and considering CO» concentration 
参数 Parameter a b c d e g h i 
re ll model (without CO;) 7.650 —0.003 0.016 —0.027 —0.036 —0.004 —0.083 —0.468 
re 12 model (with CO2) 8.022 —0.003 0.009 —0.028 —0.035 —0.004 —0.044 —0.742 
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2 未 考虑 大 气 CO 浓度 的 Irmak 模型 (rn model) 在 参数 率 定期 (a, b) 和 验证 期 (c, d) 对 冠 层 阻抗 (a, c) 和 蒸 散 量 
(b, qd) 的 模拟 效果 
Fig.2 Simulation results of canopy resistance (a, c) and evapotraspiration (b, d) with Irmak model non-considering CO; 
concentration (re i; model) in the phases of calibration (a, b) and validation (c, d) 


由 图 2 和 图 3 可 知 , 未 考虑 大 气 CO; 浓度 和 考虑 
X^ CO; 浓度 的 两 个 模型 均 能 够 较为 准确 地 反映 玉 
米 冠 层 阻抗 和 蒸 散 量 的 变化 规律 ， 模拟 值 与 实测 值 
县 有 较 高 的 一 致 性 。 在 考虑 了 LAL Ra RH, Ta ra 
U3、9 的 基础 上 , 增加 大 气 COS 浓度 影响 的 x 1 模型 
具有 更 高 的 模拟 精度 。 参数 率 定期 r. 1 模型 所 模拟 的 
冠 层 阻抗 和 蒸 散 量 与 实测 值 之 间 的 及 分 别 达 0.76 和 


0.95, RMSE 分 别 达 33.1 sm 和 34.5 Wm”, 模型 验 
证 期 冠 层 阻 抗 和 蒸 散 量 模拟 值 与 实测 值 之 间 的 R? 分 
别 达 0.68 和 0.90,RMSE 分 别 达 63.2 sm 和 49.0 W-m ^, 
模型 精度 也 较 rn 模型 略 有 提高 ,由 于 早晨 和 下 午 的 
水 汽 压 差 以 及 可 用 能 量 相 对 较 小 ， 使 得 冠 层 阻抗 变 
化 较 大 甚至 出 现 部 分 异常 值 55， 降 低 了 模型 对 冠 层 
阻抗 的 模拟 精度 ， 但 误差 均 在 可 接受 范围 内 。 
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3 考虑 大 气 CO; 浓度 的 Irmak 模型 (re n model) Zt £ 35 3E EH (a, b) 和 验证 期 (c, d) 对 冠 层 阻抗 和 蒸 散 量 的 模拟 效果 


Fig. 3 


Simulation results of canopy resistance (a, c) and evapotraspiration (b, d) with Irmak model considering CO; concentration 


(re 2 model) in the phases of calibration (a, b) and validation (c, d) 


mE, EIE TAA CO 浓度 影响 的 r np 模型 能 
够 较为 精确 地 模拟 玉米 冠 层 阻抗 对 环境 变量 的 响应 
过 程 ， 得 到 较 高 的 模拟 精度 ， 模 型 形式 简单 ， 模 型 
参数 也 较为 容易 获取 。 在 诸如 Jarvis 等 模型 中 ,包含 
许多 经 验 性 质 的 参数 ， 这 些 参数 会 随 着 植物 的 生理 
结构 .气象 要 素 的 变化 而 发 生变 化 , 但 Irmak 模型 能 
够 通过 农田 主要 微 气象 因子 较为 全 面 地 反映 冠 层 阻 
抗 对 环境 变量 的 响应 过 程 ， 从 而 实现 蒸 散发 从 叶片 
尺度 提升 至 冠 层 尺度 的 转换 ON REE r pR 
在 考虑 了 大 气 CO, 浓度 的 影响 之 后 能 够 更 为 真实 而 
全 面 地 反映 冠 层 阻 抗 对 环境 变量 的 响应 过 程 ， 实现 
尺度 转换 。 
2..0 ”模型 的 进一步 验证 

本 文 所 构建 的 两 个 模型 均 将 玉米 的 整个 生育 期 
数据 分 为 参数 率 定期 和 模型 验证 期 两 部 分 , 但 是 在 
作物 整个 生育 期 内 ,植物 的 生理 结构 、 气 象 要 素 等 
会 发 生变 化 ， 基 于 作物 生长 前 期 数据 所 建立 的 模型 
对 于 作物 生长 后 期 力 至 整个 作物 生长 季 是 否 适用 ， 
在 不 同 地 点 不 同 气象 要 素 条 件 下 ， 模 型 适应 性 如 何 ， 
模型 参数 是 否 具 有 可 移植 性 等 问题 仍 是 限制 模型 能 
否 实现 尺度 转换 和 移植 的 关键 。 为 了 检验 两 个 Irmak 
模型 的 适应 性 和 参数 可 移植 性 , 利用 通 量 观测 矩阵 
中 另外 两 个 站 点 (Site 11 和 大 满 站 )2012 年 玉米 整个 

育 期 实测 的 有 效 数据 ， 对 所 建立 的 冠 层 阻抗 模型 

做 进一步 验证 , 模拟 值 与 实测 值 对 上 比 结果 见 图 4、 
5、 图 6 和 图 7。 

从 图 4 和 图 5 可 知 ， 两 个 模型 均 能 够 较 好 地 模 
拟 Sitell 站 点 玉米 生育 期 内 冠 层 阻抗 和 蒸 散 量 对 环 


境 变 量 的 响应 过 程 ， 模 拟 值 与 实测 值 之 间 较 为 一 致 ， 
考虑 了 大 气 CO, 浓度 影响 的 r 1 模型 所 模拟 的 冠 层 
阻抗 与 蒸 散 量 模拟 值 和 实测 值 之 间 的 及 分 别 达 0.62 
和 0.93, RMSE 分 别 达 94.1 s:m FA 67.4 W:m^?, 59 
比 图 4c 和 图 5c 可 以 看 出 , x。 o 模型 模拟 值 与 实测 值 
之 间 的 误差 和 波动 性 在 整个 生育 期 内 较 小 ， 其 低估 
程度 较 reu 模型 小 , 在 6 月 中 后 期 至 8 AK, rnis 
型 所 模拟 的 ET 值 稍微 偏 低 ， 未 能 很 好 地 模拟 ET 的 
峰值 ， 相 上 比 之 下 rp 模型 很 好 地 模拟 了 整个 生育 期 
ET 的 峰值 ， 因 此 考虑 了 大 气 CO; 浓度 的 ro p 模型 能 
够 更 加 全 面 地 模拟 整个 生长 季 玉 米 的 实际 蒸 散 量 ， 
反映 其 农田 耗 水 过 程 。 
6 和 图 7 表明 未 考虑 大 气 CO 浓度 和 考虑 大 
A COo 浓度 的 两 个 Irmak 模型 均 能 够 较为 准确 地 模 
拟 大 满 站 整个 生育 期 玉米 冠 层 阻抗 和 蒸 散 量 。 未 考 
EAS CO 浓度 的 rm 模型 所 模拟 的 冠 层 阻 抗 与 蒸 
散 量 和 实测 值 之 间 的 及 分 别 达 0.57 和 0.97, RMSE 
分 别 达 45.6 sm 和 36.8 Wm”, 而 考虑 了 大 气 CO; 
浓度 影响 的 r, p 模型 所 模拟 的 冠 层 阻抗 与 蒸 散 量 和 
实测 值 之 间 的 R^ 分 别 达 0.58 和 0.97, RMSE 分 别 达 
41.5 sm ! 和 28.4 Wm”. El 6c 和 图 7c 反映 了 考虑 
XA CO 浓度 的 r, p 模型 能 够 更 加 真实 地 反映 玉米 
在 整个 生育 期 内 实际 蒸 散 量 情况 , 整个 生育 期 模型 
低估 现象 较 少 ,模拟 精度 更 高 从 6 月 6 日 到 玉米 的 
收获 期 ， 模 拟 值 与 实测 值 几乎 一 致 。 
玉米 冠 层 在 整个 生长 季 由 稀 足 向 稠密 过 渡 ， 生 
理 结构 以 及 外 界 气 象 要素 均 发 生 了 很 大 变化 ,其 蒸 
量变 化 较 大 ,尤其 是 未 来 气候 变化 情景 下 的 农田 
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图 4 REAS CO, 浓度 的 Irmak 模型 sn modeD) 对 Site 11 整个 生育 期 内 冠 层 阻 抗 (3) 和 蒸 散 量 (b, c) 的 模拟 效果 
Fig.4 Simulation results of canopy resistance (a) and evapotraspiration (ET) (b, c) with Irmak model non-considering CO; 
concentration (re 11 model) during maize growth season in Site 11 
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Fig. 5 


考虑 大 气 CO; 浓度 的 Irmak 模型 (r。 n model) 对 Site 11 整个 生育 期 内 冠 层 阻抗 和 蒸 散 量 的 模拟 效果 
Simulation results of canopy resistance (a) and evapotraspiration (ET) (b, c) with Irmak model non-considering CO; 


concentration (re n model) during maize growth season in Site 11 


蒸 散 变 得 更 为 复杂 ,因此 在 生态 系统 水 碳 平 衡 变化 
的 模拟 预测 中 ,考虑 大 气 CO 浓度 变化 对 冠 层 阻抗 
以 及 蒸 散 量 的 影响 具有 极其 重要 的 作用 。 从 图 4-7 可 
以 看 出 , 考虑 了 大 气 CO 浓度 影响 的 r o 模 型， 能 够 
更 加 精确 地 模拟 玉米 在 生理 结构 、 气象 要 素 发 生变 化 
情景 下 冠 层 阻抗 和 蒸 散 量 的 变化 ,并 且 具 有 较 好 的 
空间 可 移植 性 ,这 在 一 定 程度 上 能 够 为 种 植 结构 调 
整 、 土 地 利用 方式 和 大 气 CO, 浓度 改变 情景 下 的 玉 
米 农田 耗 水 过 程 研究 提供 一 定 的 依据 。 


2.2 ”人 敏感 性 分 析 
由 于 作物 自身 因素 以 及 外 部 环境 变量 的 改变 ， 
冠 层 阻抗 及 蒸 散 量 对 环境 变量 的 响应 也 存在 较 大 差 
F, 因此 进一步 对 改进 的 puo 模型 进行 敏感 性 分 析 ， 
引入 模型 结果 关于 环境 变量 的 敏感 性 系数 ， 如 下 式 
AR/R | OR x; 


Pm: 
| | (11) 
Axj/Xj Og R 


式 中 : s, 为 模型 结果 关于 环境 变量 * 的 敏感 系数 ， 


S 


lim 
Ax 一 0 


X; 
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Fig. 6 Simulation results of canopy resistance (r.) (a) and evapotraspiration (ET) (b, c) with Irmak model non-considering CO; 


concentration (re n model) during maize growth season in Site Daman 
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图 7 


考虑 大 气 CO; 浓度 的 Irmak 模型 (rp modeD) 对 大 满 站 整个 生育 期 内 冠 层 阻抗 和 蒸 散 量 的 模拟 效果 对 比 


Fig. 7 Simulation results of canopy resistance (r.) (a) and evapotraspiration (ET) (b, c) with Irmak model considering CO» 


无 量 纲 ， 可 以 进行 模型 不 同 环境 变量 敏感 性 的 对 
比 。 当 s, > 0 时 ， 表 示 模 型 结果 随 环 境 变量 v 的 增 大 
而 增 大 ; fis, < 0 时 ,表示 模型 结果 随 环境 变量 x; 的 


增 大 而 减 小 ; 


量 * 变化 的 敏感 程度 ， 其 值 越 大 ， 环 境 变 量 x; 
模型 结果 的 影响 越 大 。 

AR 2 和 表 3 可 知 , x. 与 R,.. RH. Us. LAL 4X 
rs 变化 方向 相反 , Rs RH Us LAL 6 r Bie 


S 


Xi 


concentration (re 2 model) during maize growth season in Site Daman 


加 而 减 小 , x 与 7, 和 大 气 CO, 浓 度 的 变化 一 致 ， 随 着 
T, REA CO; 浓度 的 增加 而 增加 。 人 敏感 系数 的 绝对 
值 大 小 表明 , x 对 Ra F RH 的 变化 最 为 敏感 ， 其 次 是 
0, re 对 Ta LAI 以 及 大 气 CO, 浓 度 的 变化 敏感 性 几 
乎 一致 ， 对 U3 和 x, 的 敏感 性 最 小 ， 表 明 在 研究 区 域 
Rn 和 RH 是 影响 六 的 主要 因子 ， 这 是 因为 在 自然 条 
件 下 ， 净 辐射 、 空 气温 度 、 湿 度 、 土 壤 水 含量 和 CO， 
浓度 是 影响 气孔 阻抗 的 主要 影响 因子 “1l。ET 与 


的 大 小 反映 了 模型 结果 随 环 境 变 
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表 2 环境 变量 变化 (+5%, 21096, 21596, 22090, 土 30%) 时 冠 层 阻力 的 变化 


Table 2 Variations of canopy resistance (re) with changes of +5%, +10%, +15%, +20%, +30% of the environmental variables 


变化 率 Change rate (96) Ra Ta RH Us ra LAI 0 CO; 
+5 -0.069 0.011 —0.056 —0.003 —0.008 —0.013 —0.040 0.012 
+10 —0.133 0.023 —0.109 —0.006 -0.017 —0.026 —0.079 0.025 
+15 -0.192 0.034 -0.158 -0.010 —0.025 —0.038 -0.115 0.038 
+20 —0.246 0.046 —0.205 —0.013 —0.033 —0.051 —0.151 0.051 
+30 -0.341 0.070 —0.288 —0.019 —0.049 —0.075 -0.218 0.079 
-5 0.076 -0.011 0.060 0.003 0.008 0.013 0.042 -0.012 
-10 0.158 —0.022 0.125 0.006 0.017 0.026 0.085 —0.024 
-15 0.248 —0.033 0.194 0.010 0.026 0.040 0.131 —0.035 
-20 0.346 —0.044 0.268 0.013 0.034 0.054 0.178 —0.047 
-30 0.569 —0.065 0.433 0.020 0.052 0.082 0.278 —0.068 


Ra: 净 辐 射 ; Ta: 气温 ; RH: 相对 湿度 ; Us 3 m 高 度 处 风速 ; rs。 空气 动力 学 阻抗 ; LAI: 叶 面 积 指数 ; 0: 土壤 有 效仿 水 量 ; COK COR 


Æo Rw net radiation; Tą: air temperature; RH: relative humidity; Us: wind speed at 3 m height; ra: aerodynamic resistance; LAI: leaf area index; 0: 


water content; CO2: air CO; concentration. 


表 3 环境 变量 变化 (+5%, +10%,+15%, +20%, 上 +30%6) 时 蒸 散 量 的 变化 


Table 3 Variations of evapotranspiration (ET) with changes of 4596, +10%, 41596, +20%, +30% of the environmental variables 


Rn Ta RH U3 Fa LAI 0 CO; 
+5% 0.054 —0.004 0.019 0.001 —0.002 0.004 0.014 —0.004 
+10% 0.110 —0.008 0.038 0.003 —0.003 0.009 0.028 —0.009 
+15% 0.166 —0.011 0.056 0.004 —0.005 0.013 0.042 —0.013 
+20% 0.222 —0.015 0.074 0.005 —0.006 0.017 0.056 —0.018 
+30% 0.337 —0.023 0.110 0.008 —0.009 0.026 0.083 —0.027 
—5% —0.054 0.004 —0.019 —0.001 0.002 —0.004 —0.014 0.004 
-1096 —0.106 0.008 —0.038 —0.003 0.004 —0.009 —0.028 0.008 
-15% —0.158 0.011 —0.058 —0.004 0.006 —0.013 —0.042 0.013 
—20% —0.208 0.015 —0.078 —0.005 0.008 —0.018 —0.057 0.017 
-30% —0.306 0.023 —0.117 —0.008 0.014 —0.026 —0.085 0.025 


R 净 辐 射 ; Ta 气温 ; RH: 相对 湿度 ; U3: 3 m 高 度 处 风速 ; rn; 空气 动力 学 阻抗 ; LAI: 叶 面 积 指数 ; 0: 土壤 有 效 含水 量 ; CO: 大 气 CO;, 浓 


Æo Rw net radiation; Tą: air temperature; RH: relative humidity; Us: wind speed at 3 m height; ra: aerodynamic resistance; LAI: leaf area index; 0: 


water content; CO»: air CO» concentration. 


Ras RH, U; LAIAU 0 的 变化 相 一 致 , ET 随 着 Ra 
RH. Us LAIFI 2 的 增加 而 增加 , ET 与 Tas ra 和 大 
A COo 浓度 的 变化 相反 ,ET 随 着 Tn r 和 大 气 CO» 
浓度 增加 而 减 小 。 敏 感 系数 绝对 值 大 小 表明 ，ET 对 
RR, 最 为 敏感 ， 其 次 是 RH 和 0, 对 7T,、LAI 和 大 气 CO; 
浓度 的 敏感 性 较 小 ， 对 U3 和 x, 的 敏感 性 最 低 , 结果 
表明 在 该 研究 区 域 R, 和 RH 是 影响 ET 的 主要 控制 
因子 ， 这 是 因为 蒸 散 发 过 程 主要 取决 于 水 汽化 所 需 
能 量 的 多 喜 ， 水 分 汽化 的 主要 动力 来 源 于 辐射 和 大 
气温 度 ， 蒸发 面 与 大 气 周 围 的 水 汽 压 差 是 决定 水 汽 
移动 的 因素 ， 而 水 汽 的 移动 过 程 很 大 程度 上 取决 于 
大 范围 的 风 以 及 蒸发 面 上 方 流动 的 空气 。 


3 结论 与 讨论 


本 研究 在 Irmak 模型 的 基础 上 , 考虑 了 大 气 
CO; 浓度 变化 对 玉米 冠 层 阻抗 的 影响 ,根据 前 人 的 


研究 结果 并 假设 大 气 CO 浓度 对 气孔 导 度 的 影响 是 
线性 关系 ,考虑 微 气象 因子 与 玉米 冠 层 阻抗 之 间 的 
关系 ， 基 于 现 有 数据 ,构建 了 考虑 CO 浓度 和 未 考 
E CO, 浓度 影响 的 冠 层 阻抗 模型 ， 并 将 其 与 Penman- 
Monteith 模型 耦合 ， 分 析 和 检验 了 所 建立 的 冠 层 阻 
抗 模型 对 环境 变量 和 大 气 CO 浓度 响应 的 模拟 结果 ， 
结果 表明 采用 LAL Ra RH, Ta rs US 以 及 土壤 
有 效仿 水量 0, 结合 大 气 CO, 浓度 能 够 更 加 全 面 地 
反映 冠 层 阻抗 和 蒸 散 量 对 外 部 环境 变量 的 响应 过 程 ， 
准确 地 模拟 玉米 在 整个 生育 期 内 半 小 时 时 间 尺 度 上 
的 耗 水 过 程 。 敏 感性 分 析 表 明 xr。 和 ET 对 R, 和 RH 
的 变化 最 为 敏感 ， 其 次 为 Ta LAI ARAS CO; 浓 
度 。 本 文 改进 的 模型 考虑 了 大 气 CO, 浓度 对 玉米 气 
孔 的 影响 ， 模 型 参数 容易 获取 ,模型 模拟 精度 优 于 
未 考虑 CO 浓度 影响 的 模型 ， 因此 该 模型 能 够 为 种 
植 结构 调整 、 土 地 利用 方式 改变 以 及 未 来 CO, 浓度 
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变化 情景 下 的 作物 耗 水 量 研究 提供 一 定 的 依据 。 

与 此 同时 ,本文 所 构建 的 模型 是 基于 “ 单 源 模 
型 ”， 对 蒸 散 量 的 模拟 具有 较 高 的 精度 ， 在 早晨 和 下 
午 由 于 水 汽 压 差 和 可 用 能 量 较 小 ， 冠 层 阻 抗 波动 幅 
度 较 大 ， 模 型 还 不 能 快速 响应 ， 使 得 部 分 冠 层 阻抗 
模拟 精度 偏 低 , 但 误差 在 可 接受 范围 内 。 其 次 , 文中 
REBT CO 浓度 上 升 对 玉米 气孔 导 度 的 影响 ， 并 
假设 CO, 浓度 对 气孔 导 度 的 影响 是 线性 的 , 但 是 并 
未 考虑 到 CO; 浓度 变化 对 玉米 叶 面 积 指数 的 影响 ， 
因此 , 需要 在 后 续 工 作 中 完善 。 此 外 ， 受 短期 观测 数 
据 限 制 , 模型 中 大 气 CO 浓度 对 冠 层 阻抗 (ro) 的 影响 
系数 需要 根据 实验 资料 做 进一步 修正 。 最 后 ， 由 于 
研究 区 域 范围 较 大 ， 虽 已 对 模型 进行 了 独立 验证 ， 
但 考虑 到 空间 复杂 性 和 异 质 性 ， 需 要 在 已 率 定 的 经 
验 模 型 的 基础 上 对 参数 可 移植 性 问题 进行 更 加 全 面 
的 考虑 。 
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